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Abstract. The fluxonium qubit has arisen as one of the most promising candidate 
devices for implementing quantum information in superconducting devices, since it is 
both insensitive to charge noise (like flux qubits) and insensitive to flux noise (like 
charge qubits). Here, we investigate the stability of the quantum information to 
quasiparticle tunneling through a Josephson junction. Microscopically, this dephasing 
is due to the dependence of the quasiparticle transmission probability on the qubit 
state. We argue that on a phenomenological level the dephasing mechanism can be 
understood as originating from heat currents, which are flowing in the device due to 
possible effective temperature gradients, and their sensitivity to the qubit state. The 
emerging dephasing time is found to be insensitive to the number of junctions with 
which the superinductance of the fluxonium qubit is realised. Furthermore, we find 
that the dephasing time increases quadratically with the shunt-inductance of the circuit 
which highlights the stability of the device to this dephasing mechanism. 
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Among the various types of superconducting qubits BEIE], the recently developed 
fluxonium qubit mu has the unique advantage of being protected against both charge 
and flux noise. This is important since both effects in general limit the performance of 
the qubits by introducing relaxation and dephasing processes. Indeed, over the last few 
years considerable effort has been made in order to understand, and consequently to 
reduce, the causes of relaxation and decoherence in different types of superconducting 
circuits 

Initially, the fluxonium qubit has been designed in order to reduce the sensitivity of 
the Cooper pair box to charge noise mu- Subsequently, it has been argued that 
it is also insensitive to flux noise [U]. In order to reach a regime in which these 
relaxation and decoherence processes are exponentially suppressed, the charging energy 
Ec = e^/2C', with the capacitance C and the elementary charge e, has to be much larger 
than the inductive scale with the inductance L and the fine-structure 

constant a = e^jhc. While it is impossible to realize such large “superinductances”, 
L 3> a~‘^C ~ lO^C, with conventional media, in the fluxonium qubit it has been 
realized by an array of (large) Josephson junctions [5]. A downside of this approach is 
that the device is then potentially plagued by spurious phase slips through the array; 
however, these have been successfully eliminated mm- Recently, different detrimental 
effects in the fluxonium qubit due to non-equilibrium quasiparticles have been addressed 
theoretically da na El ESI [IS] and experimentally im. 

Of all the processes limiting the performance of the qubits discussed above 
quasiparticle tunneling is particularly important as it is intrinsic to the superconducting 
tunneling junction and as such forms an absolute limit. While relaxation and dephasing 
mechanisms due to quasiparticle tunneling have been extensively studied in charge 
qubits [IE| EHl ESI EH E2l ES] , for the flux qubits most studies so far have concentrated 
on the relaxation due to quasiparticle processes [laiisiini. On the other hand, hrst 
investigations of the dephasing due to quasiparticle tunneling on a perturbative level in 
the junction transmission have been put forward in m El EE]. It turns out that 
treating the problem in perturbation theory leads to a diverging result due to the 
sharp peak of the quasiparticle density of states at the gap. This problem has been 
addressed by introducing as a cutoff a relaxation rate jlS] or a dephasing rate jHEE] 
that broaden the density of states, where the latter has been determined self-consistently. 
It has been discussed that in principle the divergence is lifted by treating the tunneling 
nonperturbatively as the divergence simply signals the presence of a weakly bound 
Andreev state close to the gap ig. 

Recently, it has been recognised that the presence of different nonequilibrium 
quasiparticle distributions on the different superconducting islands of the qubit (possibly 
accidentally arising during operation of the qubit and resulting in an effective 
temperature gradient of stationary nonequilibrium quasiparticle distributions) can lead 
to an additional decoherence mechanism for the flux qubits, see reference 1211 . where in 
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particular the flux qubit in the Delft-qubit design has been addressed. It has been shown 
that in this case a limitation of the dephasing time arises caused by heat currents carried 
by quasiparticles which flow through the Josephson junctions of a superconducting qubit 
as a response to the (effective) temperature gradient. 

The microscopic origin of this dephasing mechanism is the fact that a heat current 
flowing through a Josephson junction depends on the phase difference of the two 
superconductors separated by the junction. This fundamental effect has been predicted 
over 50 years ago, see pa, and has later been studied in more detail in references 
[261 1271 l28l l29] . Only very recently, the phase dependence of heat currents through 
Josephson junctions has been measured experimentally [301EIIEOI- Intriguingly, due to 
the general phase-dependence of the heat current, heat currents flowing across junctions 
of superconducting flux qubits can depend on the qubit state. This results in the 
dephasing of the qubits, since heat currents are dissipative m- 

In this paper, we investigate pure dephasing of the fluxonium qubit due to a 
non-equilibrium quasiparticle distribution. We in particular analyze the relevance of 
this intrinsic dephasing mechanism (without population relaxation) slightly away from 
the sweet spot, where it would be suppressed, see [131 HU US]- We model the non- 
equilibrium distribution of quasiparticles by an effective temperature which is different 
on different superconducting segments of the qubit. We study in details the effect of 
the heat current flowing both through the Josephson junction of the so-called black- 
sheep junction with Josephson energy T^j, which is the main source of nonlinearity, as 
well as through the array of larger junctions constituting the superinductance, which is 
shunting the black sheep. In particular, we are interested in the linear response regime 
where the difference in temperatures is small. We show that the sensitivity of the 
fluxonium qubit to heat transport and the resulting pure dephasing is suppressed by 
a factor Ef/Ej (up to logarithmic corrections) rendering the fluxonium qubit rather 
insensitive to this dephasing source. We point out that even if the effective temperature 
gradient is vanishingly small, the study of the heat conductance (rather than of the heat 
currents) is relevant: intriguingly, it can be used as a phenomenological approach for 
the investigation of dephasing due to nonequilibrium quasiparticles||] Different from our 
phenomenological approach, prior work on dephasing in flux-based qubits [131 HU HSl US] 
was based on a microscopic model. Our approach has the advantage that we are able 
to consider very small temperature gradients - corresponding to superconducting qubit 
segments with identical superconducting gaps. In this regime, the dominating effect 
of the heat-current sensitivity to the qubit state stems from the phase-dependence of a 
weak bound state originating from Andreev reflection. Due to the fact that our approach 
is nonperturbative in the tunnelling coupling, we are not plagued by divergencies as the 

f Note that a study of the charge current (or conductance) would not be appropriate here, since the 
charge of quasiparticles depends on their composition of electron- and hole-like states. In particular, this 
has as a consequence that the charge current is governed by the group velocity of electrons and holes, 
which vanishes at the gap, see e.g. |33| . This hinders the access to the phase-dependent transmission 
close to the gap, which is important for quasiparticle dephasing, via the charge current. 
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weak Andreev bound state acts as a natural phase-dependent cut-off. Additionally, our 
approach provides a nice link of this intrinsic dephasing mechanism to thermal transport 
quantities and their ability to “measure” a qubit state. 

The paper is organized as follows: In section we introduce the generic 
Hamiltonian of the fluxonium qubit. We recall the results for heat currents flowing 
through a single Josephson junctions for very small temperature gradients in section 
These findings are then used in section to investigate on the sensitivity of the heat 
current on the fluxonium states. In section]^ we show in how far the resulting dephasing 
can limit the operation of the fluxonium qubit. 

2. Fluxonium qubit 

The fluxonium is a superconducting qubit which can be thought of as a Cooper-pair 
box inductively shunted with a superinductance |1]. Its electrical circuit is shown in 
figure [^a): it consists of a Josephson junction, which we refer to as the black sheep, 
with Josephson energy Fj and charging energy Eq- This Josephson junction is shunted 
by an array of M larger Josephson junctions (with Josephson coupling energy Fj//3, 
< 1). When operated at microwave frequencies well below its self-resonant frequency 
\fEf^EcJ^/h, this array emulates a ‘superinductance’ for sufficiently large M. It is the 
presence of this superinductance which renders the fluxonium insensitive to charge noise 
while yielding a highly a n h armonic spectrum, both important points for the realisation 
of a qubit. 

2.1. Hamiltonian 

The Hamiltonian of the fluxonium consists of different parts, stemming from the 
charging and the Josephson energy of the different junctions [HTl 1^ . Charging 
effects on the black sheep junction lead to the charging energy T = AEcnf^ with 
hb the number of Cooper pairs on the capacitance plate of the black sheep. In the 
Hamiltonian, the charging energy of the black sheep serves as the kinetic energy, while 
the sum of the different Josephson energies takes the role of the potential energy. In 
particular, the potential energy of the black sheep is given by U\^ = i?j(l — cos(pb) 
with (pb the superconducting phase difference across the black-sheep junction which 
is the conjugate variable to hb = —id/difTo. For simplicity, we assume the junctions 
in the array constituting the superinductance to be equal with a somewhat larger 
Josephson energy Ej/f3, {3 < 1. The total potential energy is thus given by 17 = 
17b + M{Ej/(3){1 — cos ipm) with (pM the phase difference across each of the junctions 
constituting the superinductance. We neglect the capacitive energies of the M array 
junctions since their larger area translates into smaller charging effects. The radius r of 
the superconducting loop is supposed to be so small that the inductive energy per 
flux quantum 4>o = /ic/2e = 'neja dominates the Josephson energies i?j. In this regime, 
the magnetic flux in the superconducting loop is quantized which leads to the condition 
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Figure 1. (a) Sketch of a fluxonium qubit, consisting of a superconducting loop 

interrupted by M + 1 Josephson junctions, which is threaded by a magnetic flux 
$ « ‘ho/2- The black sheep junction has a Josephson energy Ej and a phase difference 
(/3b; the remaining M junctions with a larger Josephson energy Ej/f3 (/3 < 1) act as a 
superinductance. We consider the situation where the two electrodes separated by the 
black sheep are biased by a small temperature difference, (b) Potential energy of the 
fluxonium as a function of the phase ip for E^/Ej = 6 x 10“^ at the sweet spot, f = 1I2 
[5]. The wave functions of the three lowest lying eigenstates {Ec/Ej = 3 x 10“^) 
are depicted by red, orange and green lines, where the vertical offset indicates the 
corresponding eigenenergy. 


(/9b + MifM + 27r/ = 0 (mod27r) with the dimensionless parameter / = <h/<ho due to 
the magnetic flux $ accounting for the external magnetic held penetrating the loop. 
Using this relation, we can express the Hamiltonian in terms of one dynamical variable, 
only, which we choose to be the total phase difference across the superinductance, 
(/9 = MifM = — 27r/. Note that at fixed (/9, the phase difference (pM across each 

junction in the array constituting the super inductance becomes small for large M. As a 
result of the small phase difference across each of the large junctions [§} we are allowed 
to expand the cosine in the potential energy to second order in ipM with the result 


U = -Ej cos (/9b + -ElP^; 


( 1 ) 


(up to an irrelevant constant shift). Here, we have introduced the inductive energy of 
the superinductance 


Ej ($o/27r)^ 

/3M Leff 


( 2 ) 


with the effective inductance Leg, which increases linearly with M. It is the inductive 
term, which we obtained from expanding the cosine potential of the array to second 
order in ipM, which breaks the 27t periodicity of the total potential f/ as a function of 
(/ 9 , thus rendering the device charge insensitive mm- 


§ We here neglect potential phase slips across the array, which is justified by the large array capacitance. 
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El/E, 

Ec/E, 

Fluxonium 

6 X 10-2 

3 X 10-1 

Delft qubit 

5 X 10-^ 

2 X 10-2 


Table 1. Characteristic energy ratios for the fluxonium and Delft qubit [S5] 
when described by the Hamiltonian ([^ for the dynamic degree of freedom. The ratio 
Ej IEc approximately corresponds to the number of levels in a potential minimum, 
whereas Ej/E^ is a measure of the number of mimima of the potential which typically 
contribute to the qubit states [Hi] . 


In conclusion, the total Hamiltonian H = T + U of the fluxonium is given by 

H = AEcfi^ + - Ej cos((^ + 2nf) (3) 

with n = —id/dip = —hb the operator conjugate to p. The Hamiltonian (|^ for 
/ 1 /2 may serve as a model-Hamiltonian for different qubit types |31] • In the present 

manuscript, we mainly discuss the fluxonium qubit which is realized in the regime 
Ej S> Ec 3> El, where the latter inequality is made possible by the presence of the 
superinductance. We compare the result to the case of the Delft qubit realized in the 
regime El — Ej ^ Ec] see table for a discussion of the different energy scales. 

2.2. Qubit states 

Close to / = 1/2, the potential landscape is given by a double-well potential, which is 
an ideal starting ground to encode a qubit, see figure [^b). For the fluxonium qubit, 
we And numerically that the potential landscape features two well-localized minima for 
\^f\ ~ 0-3 with 6f = f—1/2. The two minima are situated at Pl/r- The two lowest-lying 
states in these minima, representing the qubit states, correspond to a current flowing 
clockwise/counter-clockwise in the device. However, due to the large inductance the 
corresponding currents are rather small, even though the two states differ by a large 
magnetic flux. 

The dephasing mechanism which we investigate in the following arises due to the 
fact that the different semiclassical states correspond to different superconducting phases 
on the islands which in turn makes the phase-dependent heat currents through the 
junctions dependent on the qubit state. The resulting decoherence thus projects the 
qubit on the semiclassical states Pl/r- 

We now exactly look at these semiclassical states. The position of the minima Pl/r 
is given by the solution p to the equation 
dU 

0 = — = ElP -b .Fj sin((/? -b 27r/). (4) 

dp 

At 5f = 0, the solutions are situated symmetrically around p = f) with (^9^ = —pL = p* ■ 
In general, this transcendental equation can only be solved numerically. Therefore, in the 
following, we present all analytical results in the regime El i?j, whereas the numerical 
results presented in the figures are obtained taking into account the exact solution to 
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equation|^ In the fluxonium limit El -C Ej, we have cp* ~ 7r{l—E l/E j). For small 5f 7 ^ 
0 the minima shift slightly to the right and are given by 92 ^ ^ ip* — 27r(l — tiEl/E; f)5f 
and (pL ~ - 95 * - 27r(l - 7rEL/Ej)Sf. 

In the vicinity of a local minimum cp = ^Pl/r of fhe potential shown in figure [^b), 
the potential energy can be approximated by a harmonic oscillator potential, U{ip) ~ 
U{(p) + ^Ej{ip — (p)"^ with the lowest-lying states corresponding to the ground state 
wave-functions in each of the potential minima. The spread 6(p of the wavefunction is 
given by 6(p ~ Ec /. As both the fluxonium qubit and the Delft qubit are in the 
semiclassical limit, Ec ^ i?j, we neglect the finite extent of the wave functions in the 
following and assume that they are well localized at the single value Pl/r of the phase 
variable. With that our results become independent of Ec- 


3. Heat current in a Josephson junction — linear response regime 


In the following, we want to study the effect of a small (effective) temperature gradient on 
the quantum information encoded in the fluxonium. To this end, we need to calculate the 
resulting heat current flowing through the junctions interrupting the superconducting 
loop of the fluxonium qubit. We are interested in the effect of accidental (effective) 
temperature differences; this means that we need to calculate heat currents in response 
to small temperature differences 5T, that is k-^hT/ IS. ~ 10“^. We hence evaluate the 
heat current in the linear response regime. In this section, we briefly recall the results 
of [28l |29] for heat currents in a single Josephson junction in the linear-response regime 
for small temperature gradients before we continue by investigating the effects of the 
heat current on the fluxonium qubit in the next sections. 

We consider two superconducting reservoirs with gaps Ai = A(Ti) and A 2 = A(T 2 ) 
and phase difference cj), interrupted by a junction with transmission probability D. 
Taking the number of channels contributing to transport to be N, the normal-state 
resistance of the junction is given hj R = h/{2e^ND). We now assume that the two 
reservoirs are kept at temperatures Ti = T and T 2 = T -|- ST. Being interested in small 
temperature gradients, ST/T <C 1, only, we have Ai = A 2 = A. [jj] In this regime, the 
heat current between the two superconducting electrodes is given by the linear response 
result [2HI EH] 

g(0,T,JT) = -«:(0,T)JT (5) 


with the thermal conductance 


K = 


2emkuT‘^ 


du 


ui 


Ul 


A2 


cosh^ {uj/2kjiT) (a ;2 - uj/f _ 


(o;^ — A^ cos 0 ) — sin^ — 


( 6 ) 


Here, we have introduced the energy ojc = A[1 — Dsin^( 0 / 2 )]^A of the weakly bound 
Andreev state emerging in the junction [HT]. Most importantly, the heat current depends 


II Importantly, in contrast to the procedure employed here, perturbative approaches in the tunnel 
coupling such as used in references [Ml ES] , need to introduce artificial cutoff energies to avoid emerging 
divergencies for Ai « A 2 . 
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on the phase difference (j) across the jnnction. We are interested in the tnnnelling regime 
with D 1, snch that the last term in ([^ can be neglected with the resnlt 


K = 


2e‘^RkBT‘^ 


duj 


u 


u 


A' 


cosh^ {uj/2kBT) (a;2 - 


COS0). (7) 


In the following, it is important to make the phase dependence of eqnation explicit. 
As ontlined in [Appendix A we can bring the expression in the form 


K = Ko — sin^ — In ( sin^ —]+ n ,2 sin^ — 


( 8 ) 


with the coefficients 
1 


Kq = 


2emkBT^ 

K2 = Ki|lnZi)| + 
A2 


du 


Ul 


A3 


cosh^(a;/2fcBT) ’ 


Ki = 


2e^RkBT‘^ cosh2(A/2A;BT)’ 


A 2 


UJ 


e^RkBT'^ dVaA cosh^(a;/2/cBT)(a;2 — A^) 
/•VSA 


(9) 


u 


e^RkBT^ 


OJ 


A 


(^2 _ a2 \cos\r{u;/2kBT) cosh^(A/2fcB7') 


which are independent of the phase difference. Note that, dne to the specific derivation 
that we chose for these parameters (see Appendix A), differences in the fnnctional form 

For small 


of these coefficients occnr with respect to the ones given in reference 
temperatnres compared to the critical one T <C Tc, we observe that ko,i ,2 oc 
i.e., they are exponentially snppressed at low temperatnres. 


4. Heat currents in the fluxonium qubit 

We now proceed to investigate how the phase sensitivity of the heat cnrrent through a 
Josephson junction manifests itself in the fluxonium. The situation we have in mind is 
some stationary nonequilibrium quasiparticle distribution on the fluxonium qubit. As 
mentioned above, we model this distribution by an effective temperature, in order to 
keep the discussion simple. [J] In particular, we treat the situation where one side of the 
black-sheep Josephson junction is at an elevated temperature T 2 = T + ST with respect 
to the other which is at temperature Ti = T, see figure [^a). As a result quasiparticles 
tunnel from the “hot” to the cold reservoir. As seen in the last section, the resulting 
heat current depends on the superconducting phase differences and thus on the state of 
the qubit. As we have assumed the M Josephson junctions of the array emulating the 
superinductance to be equivalent, the temperature gradient is distributed among the 
M elements with the temperature difference STm = ST/M on a single junction. The 
heat current Q flowing into the cold reservoir, which is held at temperature Ti = T (or 
equivalently in the heat current flowing out of the hot reservoir, kept at temperature 


^ Note that more general nonequilibrium quasiparticle distributions can be analyzed by replacing the 
Fermi functions (entering through the cosh term in equations^to^ by arbitrary distribution functions. 

















9 


Dephasing of fluxonium qubit 

T 2 = T + (5T Q is given by the sum of two terms 

Q = T, 6T) + QiipM, T, ST/M), (10) 

where the hrst term is the contribution of the black sheep and the second term is due 
to the large Josephson junction connecting the array of the superinductance to the cold 
reservoir. 

The decoherence of the fluxonium qubit is triggered by the difference 5Q = 
IQr ~ Qi] of the heat currents flowing in the device when the qubit is in the state 
R/L with if = Following our previous work, [23], we introduce the sensitivity 

^ 6Q 
IQr + Ql\ 

as a measure of correlation between the heat current and the qubit state. As a large 
sensitivity corresponds to a large difference of heat currents for different qubit states, 
we expect that this in turn leads to fast dephasing of the qubit; an expectation which 
we conhrm below. However, in a hrst step, we want to calculate the sensitivity for the 
huxonium. 



4-1. Effect of the number of junctions implementing the superinductance 


In this section, we show that the sensitivity of the heat currents to the huxonium 
state only depends on the ehective parameters El and Ej of the Hamiltonian, given in 
equation]^ and not on the specihc number M of junctions (or their asymmetry factor (3) 
with which the superinductance Les is emulated. In order to realise a huxonium qubit 
modelled by the Hamiltonian given in equation ([^, a large amount M of array junctions 
is needed, which all have a Josephson energy larger than the one of the black sheep by a 
factor In view of equation ([^, to obtain a given value of El, the coupling strength 
of each array junction needs to scale like M oc /3~^. 

It is clear that from the two terms in equation (10), only the one stemming from 
the junction of the array forming the superinductance, Q{(Pm, T, 6T/M), could possibly 
depend on M. The variables occurring in the argument of the heat current depend on 
M via (fM = —(v5fe + 2vr/)/M and 5T/M. Exploiting equations ([^ and ([^, we hnd that 
for large M, where (fM 0, the only relevant term for Q{fM ~ 0,T,ST/M) is given 
by the phase-independent part of the thermal conductance of the array junction, which 
we denote by Kom(T). As a consequence, the state-dependent heat current difference 
6Q = IQr — Ql\ depends only on the heat current through the black sheep. Hence in 
the following only the dependence on M of the term kqm{T) entering the denominator 
of the sensitivity has to be investigated. 

The thermal conductance Rqm is proportional to Rfj, the normal-state resistance 
of the outer junction of the array. Now, due to the generalised Ambegaokar-Baratoff 
relations [38], the normal-state resistance is inversely proportional to the Josephson 


Being interested in the stationary situation, we can make this assumption as long as coupling to 
phonons is neglected. 
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Figure 2. Sensitivity as a function of El/E} for D — 10“^ and k^T = 0.1 A. 
The full result of equation (11) (full blue line) is compared to its approximation for 
SfjE^/E} —0 of equation (12) (red dashed line). 


energy of the corresponding junction and we thus hnd kqm = with kq the phase- 

independent contribution to the thermal conductance of the black sheep. Therefore, 
Q{(pM ~ 0,T,6T/M) = —f3~^M~^Ko6T = —{El/Ej)koST is independent of M, due 
to the cancellation of the factors M and [3 occuring in kqm and 6Tm. Moreover, as 
El <C Ej, the heat current through the superinductance is negligible compared to the 
one through the black sheep, which is proportional to —kST. We therefore completely 
neglect the heat current through the superinductance in the following. 


4-2. Sensitivity of the heat current to the state of the fluxonium qubit 


Having found that the sensitivity is independent of the specihc realization of the 
superinductance, we here present an analytical expression for the sensitivity s in the 
regime 5/ -C 1 and El Ej, which is relevant for the fluxonium qubit. We have seen 
in section 2.2 that ifR = ip* — 27r(l — tiEl/ EflSf and pl = —p* — 27r(l — ttEl/ EflSf 


with p* = 7r(l — El/Ej). Evaluating s of equation [ll] to Erst order in 6f and leading 
order in El/E j yields the hnal result 

27r^Ki 


Ko 




( 12 ) 


This expression shows that the sensitivity depends quadratically on the ratio El/E j with 
logarithmic corrections. Thus, the heat current in the fluxonium qubit is less sensitive 
to the qubit state than it is the case for the Delft qubit, due to the smaller ratio El/E}. 
This comparison can be easily performed, since a description of both qubits is possible 
following equation (|^, with the effective qubit parameters given in table In what 
follows, we show that the fluxonium qubit therefore enjoys an increased protection with 
respect to quasiparticle processes as compared to the Delft qubit. In hgure we have 
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plotted the approximate result of the sensitivity of equation exhibiting a quadratic 
dependence for small E^/Ej; A comparison to the exact expression [IT| indicates that 
the approximation is valid for El/Ej < 0.05. 


5. Dephasing time 

In a recent work m, we have demonstrated that the sensitivity of the heat current in 
a flux qubit (in the Delft qubit design) leads to a dephasing of the qubit. The reason 
for this dephasing is the fact that for non-zero sensitivities the tunnelling probabilities 
of quasiparticles depends on the state of the qubit; hence quasiparticles which tunnel 
through the Josephson junction, can dephase the qubit. Importantly, the heat current 
incorporates both the phase-dependent quasiparticle transmission probabilities through 
the junction as well as the quasiparticle distribution functions. This forms the basis 
for our argument that the heat current captures the relevant properties leading to 
qubit dephasing due to quasiparticle tunneling. We therefore propose to investigate this 
transport property, in order to access quasiparticle dephasing on a phenomenological 
level. We gain additional conhdence in our results for the qubit dephasing obtained from 
this phenomenological approach by noticing that in the regime of large temperatures 
(where the cut-off due to the Andreev bound state becomes unimportant) our approach 
reproduces the perturbative results from the microscopic model 


Following reference m, see also [Appendix B| we can derive the expression 


^ = 


fELSf 


doo 


ui 




V J Ja + A"^) cosh‘^{uj/2kBT) 

for the inverse dephasing time, valid to lowest order El/Ej and 5/; here, u* = 
A(1 — D sin^ ^)^/^ is the bound state energy at the phase difference (p* at 6f = 0. Note 
that at the sweet spot the heat currents are equivalent for the two qubit states which 
results in the vanishing of the dephasing rate, up to exponentially small corrections in 
Ej/Eq, which are not considered in this paper. 

Owing to our phenomenological approach, the dephasing time can be directly 
brought into contact with the sensitivity of the heat currents flowing in the device 
to the qubit state. In order to estimate this link, we consider the low temperature 
regime, T,6T ^ A/ks, and write down the product between dephasing time and the 
difference in heat currents in the two qubit states, 6Q. This function gives us an idea 
about the energy which is transferred by the difference of heat currents in the two qubit 
states in the time, which the qubit needs to dephase. It turns out that to lowest order 
in El/Ej and 6f, we obtain the simple result 

^ ^ 4:kBT^f^ ^ ’ 

The full expressions are presented in Appendix C Equation 14 shows that the value of 
t^6Q only depends on the detuning 5f from the sweet spot and is otherwise independent 
of any of the qubit parameters. The additional parameters A, T, 5T entering the 
expression describe the heat current due to the flow of the quasiparticles. 
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Figure 3. Dephasing time, r^, multiplied by the level splitting of the fluxonium e and 
by 6f as a function of temperature (full blue line). The transmission of the junction is 
assumed to be D = 10“^, Ei^/Ej = 6 x 10“^ and the number of channels contributing 
to transport N « 10^. The red dashed line shows a fit d x exp(A/fcBT) with d = 10“^. 


Starting from equation 14, we only need to find 6Q in order to obtain the dephasing 


time that limits the performance of the superconducting qubit. In the limit 6f small, 
the two states I-^l) and I-^r) have a similar heat current with Ql + Qr = 2Q{ip = ip*), 
such that we obtain the expression 6Q = 2s\Q{ip = p>*)\- Additionally, in the limit 
El <C £^j, p>* is close to vr and thus = tp*) = —n^bT which leads to the final result 

1 

s 

We see that as expected the dephasing time is inversely proportional to the sensitivity 
of the heat current that is flowing through the structure to the state of the qubit. As 
we have seen before, s is proportional to (El/Ej)'^ (up to logarithmic corrections); 
consequently the fluxonium qubit has the beneht of a larger dephasing time than the 
Delft qubit. Different from the dephasing caused by charge and flux noise, the inverse 
dephasing time due to the transport of quasiparticles is not exponentially suppressed 
in the fluxonium regime El Ej. However, due to the fact that Kq oc exp^—A/ksT), 
an exponential improvement of the dephasing time limited by quasiparticle tunneling 
can be reached by lowering the temperature. This is expected since it corresponds to 
decreasing the quasiparticle occupation. 

In order to make a quantitive statement relevant for applications, it is useful to 
introduce the dimensionless number er^/h with e: the level splitting of the qubit, rather 
than looking at the dephasing time only. In particular, since the inverse level splitting, 
yields a measure of the time of a qubit operation, the parameter er^/h indicates the 
typical number of single qubit gates which can be performed before coherence is lost. We 
numerically determine the level splitting from the difference of the energies of the ground 
and excited state from the Hamiltonian given in equation]^ In figurewe show a plot of 
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the function 5f /h. Since with small changes in /, the level splitting increases linearly 
with 6f, while is proportional to we choose to multiply the parameter of interest 
with Sf in order to get a function which is essentially independent of Sf. Figure [^shows 
a logarithmic plot of er^Sf/h for small values of the temperature, T <C A//cb- The 
ratio between the two time-scales, and h/e, occurs to be approximately exponentially 
suppressed with increasing temperature. This can be seen by using equations 12 and 
15, where we see that oc Moreover, from the dehnitions of Ki{T) given 


in equation]^ for small temperatures we have T‘^ki{T) oc exp(— A/ZcbT). 

Finally, hgure shows that for small values of the temperature it is possible to 
perform a great number of operations on the qubit state before they become unreliable 
due to dephasing of the two-level system. For example for 6f = 10“^ and k-^T < 0.15A, 
we have that er^/h > 10^. 


6. Conclusion 

In this paper, we have studied the impact of the phase sensitivity of the quasiparticle 
transport on the coherence properties of the fluxonium qubit. Using a phenomenological 
approach, based on the study of heat currents carried by quasiparticles and the 
associated heat conductance, we have shown that the dephasing time is inversely 
proportional to the sensitivity, a quantity describing to which extent possible heat 
currents flowing in the device depend on the state of the qubit. We have shown that 
the sensitivity of the heat current to the qubit state depends quadratically on the ratio 
of the characteristic energies of the fluxonium qubit but not on the number of 
junctions with which the superinductance, a relevant ingredient of the fluxonium qubit, 
is realised. The independence of the number of array-junctions can qualitatively be 
traced back to the fact that, at small temperature gradients, the heat current in the 
arms of the loop constituting the fluxonium qubit is mainly given by the heat current 
in the so-called black sheep junction, which does not depend on M. The fact that the 
sensitivity can be reduced by lowering the ratio El/Ej has the important result that 
the fluxonium qubit is less affected by dephasing due to quasiparticle tunnelling through 
the Josephson junction when compared to the Delft qubit design. We furthermore hnd 
that the dephasing mechanism is exponentially suppressed with temperature due to 
its origin in quasiparticle tunneling. However, we have shown that at moderately low 
temperatures the resulting dephasing time is demonstrated to be large enough to easily 
allow an excess of 10^ operations before the qubit dephases. 
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Appendix A. Evaluation of the linear-response coefficients for weak tunnel 
coupling 

We aim to get an understanding of the phase dependence of the linear response coefficient 
K of the heat current in the weak tunnel coupling regime, H -C 1 as given in equation 


K = 


duj 


u 


u 


A2 


-(o;^ — cos 0), (A.l) 


2emk^T^ cosffi {u/2kBT) - a;2)2 

by rewriting it in terms of the expression given in equation In this appendix, we 
outline the procedure, which we apply to obtain the coefficients given in equation 

As a starting point, it is useful to evaluate the logarithmic divergence occuring in 


the linear response coefficient shown in equation A.l Since the divergence stems from 
values of u in the vicinity of the superconducting gap A, this logarithmic divergence can 
conveniently be extracted by setting ca ~ A in all contributions of equation A.l with a 
smooth dependence on u in the vicinity of A. This consideration leads us to determine 
the integral 


A3 sin^ t 


r-VSA 


^div — 


e^Rk^T^ cosh" (A/2fcBT) 


duj 


ui 


A3 sin" I 


2e^RkBT^ cosh" {A/2kBT) 


In 


a;2 — A2 (l — H sm 

n • 2 

U sm — 


2 ± 


1 ) 


(A.2) 


The latter is the only divergent contribution in equation A.l With R~^ oc D, we find 
that the leading contribution to this term is of the order D\n.D. The remaining part of 
the integral can safely be expanded for small Hsin" We find the remaining part by 
simply substracting the logarithmic divergence from the full coefficient Krem = k — K^n, 
leading to 




2emkBT^ 


du 


u 




-A" 


cosh" {pj/2kBT) (a;2 — 


(ca" — A" COS0) (A.3) 


A3 sin" t 


r>x/3A 


emkBT^ cosh" (A/2fcBT) 


du 


u 


a;2 — A2 (l — H sin" |) 


Indeed, the expansion of equation A.3 to linear order in D shows no divergent behavior 
anymore. We find 




2e‘^RkBT 

A" 


du 


u 


cosh" {u/2kBT) 

~ , a;" 

du- 


sin"f 


2e^RkBT cosh" {u/2kBT) u"^ — A^ 


A" 


-Vs A 


du- 


u 


u 


A 


. 2 

sm — 
2 


2e‘^RkBT a;2 — A2 \_cos\Y {u/2kBT) cosh" (A/2 A;b7"). 

This expression, together with the contribution from the logarithmic divergence Kdiv, 
yields the results presented in equations [^andj^in the main text. 
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In order to investigate the impact of possible temperatnre gradients on the dephasing 
of the fluxonium qubit, we follow the lines of reference m- We therefore use a model 
Hamiltonian, H^od = Ho + Hi, where Hq describes the qubit as a two-level system with 
states IV'l) and |V’r) and normal quasi-particle reservoirs at different temperatures. The 
coupling between them is given by Hi. More specifically, we have 


hfmod 


2^ 


/=1,2 A:,£7 




ka^l,ka 


I + H.c. 


(B.l) 


k,q,a 


The matrices D,j = 0,3 are Pauli matrices in the qubit space. The level splitting 
between the qubit states is given by £; coupling between them is supposed to be weak 
and is neglected here. In the reservoirs, / (/ = 1, 2), the creation (annihilation) operators 
of particles with momentum k and spin a are given by It is the state- 

dependent coupling between qubit and reservoirs occurring in the interaction part of the 
Hamiltonian, V 0/3 = (Va ± ^l)/ 2, together with the density of states of the reservoirs, 
which takes account for the phase dependence (and hence for the dependence on the 
fluxonium states) of the heat current due to the superconductors. 

We recover the state-dependent qubit-reservoir coupling and the density of states by 
comparing the heat current obtained from the model Hamiltonian in the linear response 
regime 


/~)mod _ 

VR/L — 


n5T 


AhkeT^ 


du 


ui 


cosh^(a;/2fcBT) 


R/L R/L 


Hn 


(B.2) 


to the one calculated starting from equations (|^ and (10) for the black sheep. Note 
that we here neglect the effect of the heat current through the array junctions, since 
no relevant dependence on the qubit state occurs there. In order to extract the local 
density of states of the continuum states above the gap, we use the relation 

|a;|(a;2-A2)V2 






■Iv^r/l) 


= n?- 


|a;|(a;2-A2)V2 


U‘ 


— A 2 (l — D sin 


2 V^R/L ' 


Here, is the normal conducting density of states including spin. Using this form for 
the density of states we have the following expression for the tunnelling matrix elements 
in the tunnelling regime, D -C 1 , 



^12 



A2 

- cos (fn/L 


CO 


(B.3) 


Here, V 12 is the tunnelling amplitude of the junction figuring as the black sheep in the 
fluxonium qubit. It is linked to the normal state resistance by i? = ^/(Tre^n'j’n^Rl). 

With the help of this simplified model, we proceed to study the dynamics of the 
qubit state. Starting from the density matrix of the full system consisting of the qubit 
coupled to reservoirs, we trace out the reservoir degrees of freedom and write down a 
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master equation for the reduced density matrix of the qubit, p{t). It is helpful to rewrite 
the interaction part of the Hamiltonian as Hi = Pr-Br + PlBi, with the projectors on 
the qubit states, a = R, L, 

Pa = \‘lpa){i’a\, Ba = Va^ + H.C. (B.4) 

k,q,a 

Following a standard procedure, see for example Ref. [M], the master equation for the 
density matrix of the qubit then takes the form 

p{t) = -^[Hs,pit)]+ 7a/3(^/3P(t)Pa - ^{PaP/3,P(t)}) 

a,/3=R,L 


with the transition rates between qubit states a, {3 = R, L, 


= 2 / {{B^(t),Bi,(0)})dt. 


(B.5) 


The relaxation behavior of the qubit becomes particularly clear when rewriting the 
master equation in terms of a Pauli rate equation for the pseudo-spin states of the 
qubit, S{t) = Tr[p(f)T] = [plrW + Prl(^),*PLR( t) - Wl(^),Pll(^) - Prr(^)]'^- We 
obtain from equation (B.5), 

S{t) = S{t) xh- 0)T (B.6) 

The Pauli rate equation ( B.6| ) contains a precession of the pseudospin around a pseudo- 
magnetic field, h = {0,0,e/h)^, determined by the level splitting between qubit states. 
Most importantly, there is also a relaxation of the coherences of the reduced density 
matrix with the dephasing rate given by 


.-1 


Ti 




27rl + 7ll)- 
It is found to have the explicit form. 


.-1 


TT 


Tj, 


= — duj 
2h ./a 


I4(nfn, 


2e2p 


du 


R„R^l/2 




- VL{nin 

- A2 


WLp/2 


cosh^(a;/2/cBP) 


(B.7) 

cosh“^(ci;/2fcBP) (B-8) 

cos<pl — \/a;^ — A^ cos <pr j 


where we introduced cn* = A{1 — D sin^ 


Appendix C. Link between heat currents and dephasing time 


As discussed in the main text, the dephasing time can be directly brought into 
connection with the heat current flowing through the qubit due to a finite temperature 
gradient and its sensitivity to the qubit state. In order to estimate this link, we consider 
the linear response regime 5T <^T and write down the product between dephasing time 
and the difference in heat currents in the two qubit states, 5Q, 

cos(Pr| cosh“^ 


t^5Q = 


A^6T da;a;2^J^J^|cos (^l 


2)2 


2kBT 


2fcBT2j“ duj ( - A^ cos - a/cj^ - A^cos^^r^ cosh ^ 


-.(C.l) 


2kBT 
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As we noticed before, the main contribution to these integrals stems from contributions 
of u) close to the superconducting gap A. We hence introduce a; ~ A in those factors 
which are smooth functions of u in the vicinity of A 

cos (Pr| cosh“^ ^ 




'<t>SQ = 


A^6T dw A2 ^J^^|cos(Pl 


2kBT 


2ksT‘l 


°° rh . 


(y\/A? — A? cos</9l — i/A^ — A^ cos </9 r j cosh ^ 


-(C.2) 


2kBT 


and observe that this leads to a cancellation of the integral terms in numerator and 
denominator of equation C.2[ This underlines the close connection between the heat 
current sensitivity to the qubit state and the occurring dephasing mechanism. Further 
simplifying we find 

A^5T I cos (pL — cos (pR I 


t^5Q = 


Ak^T'^ (I sin %| — | sin 


(C.3) 


In lowest order in El/Ej and Sf, this leads to the result presented in equation 14 
the main text. 
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